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SENSOR ARRAY INDUCED PHASE ANGLE CALCULATION BASED ON MIXED 

SIGNALS 

TECHNICAL FIELD 

The invention relates generally to signal processing and more particularly to 
5 demodulation of signals from fiber optic sensor arrays. 

BACKGROUND 

Fiber optic sensor arrays of a time division multiplexed ("TDM") system are often 
used to measure a change in an external parameter, for example, acoustic vibration, fluid 
pressure variations, acceleration, and magnetic field intensity. The fiber optic sensor array 
10 employs a phase generated carrier with a period T to measure the change in the external 
parameter at a given sampling rate. The fiber optic sensor array converts a time-varying 
phase angle associated with the external parameter to an amplitude variation on an output 
pulse of light. 

The phase angle is measured through various demodulation techniques of the output 
15 pulse. Typical demodulation techniques employ a quadrature component Q and an in-phase 
component I of the output pulse. The quadrature component Q corresponds to a sine of the 
phase angle, and the in-phase component I corresponds to a cosine of the phase angle. An 
arctangent of the ratio Q/I is equal to the phase angle. The magnitude of the change in the 
external parameter can then be calculated from the change in the phase angle. 
20 Calculation of the quadrature component Q and the in-phase component I requires 

multiple samples of the output pulse at specific intervals of the phase generated carrier. The 
phase generated carrier comprises a period that is significantly longer than a period of the 
output pulse. The longer period of the phase generated carrier requires the samples to span 
several output pulses to obtain each required interval of the phase generated carrier. It is 



desirable to increase the rate of calculation of the phase angle to promote an increase in 
throughput. For a constant sampling rate, the number of required samples must be reduced to 
increase the rate of calculation of the phase angle. 

Thus, a need exists for reduced numbers of required samples for demodulation 
5 techniques of fiber optic sensor arrays that employ phase generated carriers. 

SUMMARY 

The invention in one embodiment encompasses a method. A sensor array employs a 
parameter to induce a time-varying phase angle <p on an optical signal that comprises a phase 
generated carrier based on a single frequency. The time-varying phase angle cp is calculated 

10 through employment of a plurality of samples obtained approximately at a time t s? where the 
plurality of samples are based on a plurality of mixed signals output from a mixer component 
that employs an input signal based on the optical signal of the sensory array. 

Another embodiment of the invention encompasses an apparatus. A sensor array 
employs a parameter to induce a time-varying phase angle (p on an optical signal that 

1 5 comprises a phase generated carrier based on a single frequency. The apparatus comprises a 
processor component that calculates the time-varying phase angle cp through employment of a 
plurality of samples obtained approximately at a time t s , where the plurality of samples are 
based on a plurality of mixed signals output from a mixer component that employs an input 
signal based on the optical signal of the sensory array. 

20 A further embodiment of the invention encompasses an article. The article comprises 

one or more computer-readable signal-bearing media. A sensor array employs a parameter to 
induce a time-varying phase angle cp on an optical signal that comprises a phase generated 
carrier based on a single frequency. The article includes means in the one or more media for 
calculating the time-varying phase angle (p through employment of a plurality of samples 



obtained approximately at a time tg, where the plurality of samples are based on a plurality of 
mixed signals output from a mixer component that employs an input signal based on the 
optical signal of the sensory array. 

DESCRIPTION OF THE DRAWINGS 

5 Features of exemplary implementations of the invention will become apparent from 

the description, the claims, and the accompanying drawings in which; 

FIG. 1 is a representation of one exemplary implementation of an apparatus that 
comprises one or more lasers, one or more optical switches, one or more phase modulators, 
one or more sensor arrays, one or more photodetectors, one or more mixer components, one 
10 or more oscillator components, one or more filter components, and one or more processor 
components for performing a calculation of a phase angle. 

FIG. 2 is a plot of an exemplary interference pulse, an exemplary impulse function, 
and an exemplary filtered pulse with an exemplary filter component of FIG. 1 . 

FIG. 3 is a plot of an exemplary roll-off function of the exemplary filter component of 

15 FIG. 2. 

FIG. 4 is a plot of one or more exemplary filtered signals with the exemplary filter 
component of FIG. 2. 

FIG. 5 is a plot of an exemplary in-phase term I(t, (p) and an exemplary quadrature 
term Q(t, q>) as functions of the phase angle for the filtered signals of FIG. 4. 

20 FIG. 6 is a plot of an exemplary error for the calculation of the phase angle of FIG. 5. 

FIG. 7 is a plot of an exemplary pulse height deviation for filtered signals of FIG. 5. 



DETAILED DESCRIPTION 

Turning to FIG. 1, an apparatus 100 in one example comprises a plurality of 
components such as computer software and/or hardware components. A number of such 
components can be combined or divided in the apparatus 100. An exemplary component of 
5 the apparatus 100 employs and/or comprises a set and/or series of computer instructions 
written in or implemented with any of a number of programming languages, as will be 
appreciated by those skilled in the art. 

Referring to FIG. 1, the apparatus 100 in one example comprises one or more lasers 
102, one or more optical switches 104, one or more phase modulators 106, one or more 

10 sensor arrays 108, one or more optical receivers 1 10, one or more mixer components 111, one 
or more oscillator components 112, 113, 114, and 115, one or more filter components 116, 
and one or more processor components 118. In one example, the apparatus 100 demodulates 
an optical signal to measure a change in a parameter, as described herein. The laser 1 02 in 
one example comprises a continuous wave laser. The laser 102 generates and sends an 

15 optical signal through the optical switch 104 and the phase modulator 106 to the sensor array 
108. 

The optical switch 104 in one example comprises a time division multiplexed 
("TDM") switch. The optical switch 104 gates the optical signal such that the optical signal 
comprises a stream of optical pulses. The phase modulator 106 impresses a phase generated 

20 carrier ("PGC") 120 on the stream of optical pulses. For example, the laser 102, the optical 
switch 104, and the phase modulator 106 cooperate to create one or more optical pulses 122 
that comprise the phase generated carrier 120, as will be understood by those skilled in the 
art. The optical pulse 122 comprises a period T pu i se . The period T pu i se in one example is 
approximately between 50 nanoseconds and 1000 nanoseconds. The phase generated carrier 

25 120 in one example is based on a single frequency f pgc and comprises a modulation depth of 



M. A period T pgc comprises a relationship with the frequency f pgc = 1 / T pgc , as will be 
understood by those skilled in the art. The frequency f pgc in one example is approximately 
between 20 MHz and 200 MHz. The phase generated carrier 120 is associated with a 
demodulation phase offset p. The phase generated carrier 120 creates a time-varying phase 



5 angle equal to M • sin 



T 

V ps c J 



The sensor array 108 in one example comprises one or more sensors 130, 132, and 
134, for example, mismatched path interferometers. The sensor array 108 splits the optical 
pulse 122 into one or more optical pulses 124, 126, and 128, for example, one pulse per 
sensor. The optical pulses 122, 124, 126, and 128 in one example are substantially the same. 

10 The sensors 130, 132, and 134 of the sensor array 108 receive the optical pulses 124, 126, 
and 128, respectively. The sensors 130, 132, and 134 of the sensor array 108 in one example 
employ one or more parameters and the optical pulses 124, 126, and 128 to create one or 
more respective interference pulses 136, 138, and 140. Exemplary parameters comprise 
acoustic vibration, fluid pressure variations, acceleration, and magnetic field intensity. For 

15 example, the sensor 130 splits the optical pulse 124 into a first portion and a second portion. 
The sensor 130 employs the parameter to induce a time- varying phase angle (p on the first 
portion of the optical pulse 124, relative to the second portion of the optical pulse 124. The 
sensor 130 recombines the first portion of the optical pulse 124 with the second portion of the 
optical pulse 124 to create the interference pulse 136. A time- varying amplitude variation of 

20 the interference pulse 136 represents the time-varying phase angle cp between the first portion 
and the second portion of the optical pulse 124. 

The optical pulses 122 comprise an intermediary spacing such that the interference 
pulses 136, 138, and 140 comprise a relatively small spacing, for example, a high duty cycle, 
as described herein. The interference pulses 136, 138, and 140 comprise a period 



substantially equal to the period T pu i S e of the optical pulse 122. The sensor array 108 sends 
the interference pulses 136, 138, and 140 to the optical receiver 110 in a pulse train 142, for 
example, in a serial fashion. For example, the optical pulse train 142 comprises the 
interference pulses 136, 138, and 140. 
5 The optical receiver 1 10 in one example comprises one or more photodiodes 144. In 

a further example, the optical receiver 110 comprises a transimpedance amplifier 146. The 
optical receiver 110 in one example comprises a polarization diversity receiver system as 
disclosed in U.S. Patent No. 5,852,507 to Hall, issued May 10, 1996, and assigned to Litton 
Systems Inc. The optical receiver 110 receives the optical pulse train 142. The optical 
10 receiver 110 then creates one or more respective analog electrical signals that represent the 
interference pulses 136, 138, and 140 from the optical pulse train 142. For example, the 
optical receiver 1 10 converts a magnitude of power of the optical pulse train 142 to a voltage 
signal. The analog electrical signals in one example are defined by s(t, M, 9), as described 
herein. 

15 The mixer component 1 1 1 in one example comprises a high-frequency analog mixer. 

The mixer component 1 1 1 in one example employs an input signal to create a plurality of 
mixed signals. The input signal is based on the optical signal, for example, the input signal 
comprises the analog electrical signals that represent the interference pulses 136, 138, and 
140 from the optical pulse train 142. The mixer component 111 mixes the input signal with 

20 one or more sinusoidal signals from the oscillator components 112, 113, 114, and 115 to 
create the plurality of mixed signals. The mixer component 111 then sends the plurality of 
mixed signals through the filter component 116 to create a plurality of filtered signals, as 
described herein. 

The oscillator components 112, 113, 114, and 115 in one example comprise 
25 waveform generators. The oscillator components 112, 113, 114, and 115 generate the one or 
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more sinusoidal signals for the mixer component 111. The sinusoidal signals in one example 
comprise: 

I s (t) = sin [47cf pgc (t + tb)L 

Ic(t) = COS[47lf pgc (t + to)L 

5 Q s (t) - sin [ 27if pgc ( t + to ) ], and 

Qc(t) = COS[27Tfp gc (t + to)]. 

For example, the oscillator component 112 generates I s (t), the oscillator component 113 
generates Ic(t), the oscillator component 114 generates Q s (t), and the oscillator component 
115 generates Q c (t). A phasing difference indicated by a time to between the sinusoidal 

10 signals and the input signal is arbitrary and not controlled, as will be appreciated by those 
skilled in the art. For example, the calculation of the phase angle cp is independent of the 
phase offset of the phase generated carrier with respect to the sinusoidal signals. 

The filter component 116 comprises a low-pass filter. The filter component 116 in 
one example filters the plurality of mixed signals to create a plurality of filtered signals, as 

1 5 described herein. An ideal low-pass filter performs time integration over a time period T. 
For example, an ideal low-pass filter filters the mixed signals to create the filtered signals: 

F is (M,q>)= js(/,M,0>)-/ 4 (O<ft = BTJ 2 (M) cos q> sin (47if pgc to), 
F ic (M,<p)= \s(t,M,<p)-I c (t)dt = BTJ 2 (M) cos 9 cos (47cf pgc to), 
F qs (M, q>) = js(t,M 9 <p) ■ Q s {t)dt = BTJi(M) sin cp cos (27if pgc to), and 

20 F qc (M, cp) = Js(* , M, <p) ■ Q c (t)dt = - BTJi(M) sin cp sin (2Tcf pgc to), 

where Ji(M) comprises a first order Bessel function and J 2 (M) comprises a second order 
Bessel function. The time integration comprises a time period T such that time-varying 
components of the mixed signals average out to zero, as will be understood by those skilled 
in the art. The filter component 116 in one example comprises a Gaussian filter, a fourth 



order Bessel filter, or a fourth order real pole filter that approximates an ideal low-pass filter, 
as will be appreciated by those skilled in the art. The time to is dependent on a time of travel 
for the pulse 122 from the optical switch 104 to the optical receiver 110; therefore, the 
quantity 27rf pgc to changes slowly, as will be understood by those skilled in the art. 
5 The processor component 1 18 in one example comprises a digital signal processor. In 

a further example, the processor component 1 1 8 comprises an analog-to-digital converter 
component 148. The processor component 118 in one example comprises an instance of a 
computer-readable signal-bearing media 150, as described herein. The analog-to-digital 
converter component 148 converts the filtered signals from the filter component 116 into 

10 digital signals. The processor component 1 18 in one example serves to sense a change in the 
parameters by employing the time-varying amplitude variation of the interference pulses 136, 
138, and 140 to calculate the time- varying phase angle <p. 

An illustrative description of exemplary operation of the apparatus 100 is presented, 
for explanatory purposes. The laser 102, the optical switch 104, and the phase modulator 106 

15 cooperate to create the one or more optical pulses 122. The sensor array 108 splits the optical 
pulse 122 into the optical pulses 124, 126, and 128. The sensors 130, 132, and 134 employ 
the parameters and the optical pulses 124, 126, and 128 to create the interference pulses 136, 
138, and 140. The sensor array 108 sends the interference pulses 136, 138, and 140 as the 
optical pulse train 142 to the optical receiver 110. 

20 The optical receiver 1 10 creates an analog electrical signal that represent the one or 

more interference pulses 136, 138, and 140. For example, the analog electrical signal is 
defined as s(t, M, cp): 

s(t,M, <p) = A + B- cos(M ■ sin(2^ gc ■ t)- <p) 

= A + B • cos (p • cos(M • sin(2^T ■ t)) + B • sin (p • sin(M • sin(2^T • t)) , 



where A is an average signal level, B is an interference term signal level, M is the 
modulation depth, f pgc is the frequency of the phase generated carrier, and cp is the phase 
angle. The phase angle of s(t, M, cp) comprises a first portion due to the phase generated 
carrier, M • sin(2;{f • t), and a second portion due to the parameter, cp, as will be understood 

5 by those skilled in the art. 

The mixer component 111 mixes the analog electrical signal with the sinusoidal 
signals I s (t), ^(t), Q s (t), and Q c (t) to create the mixed signals M is (t, M, cp), M ic (t, M, 9), M qs (t, 
M, cp), and Mq C (t, M, cp), respectively. The mixer component then sends the mixed signals 
Mj s (t, M, cp), Mi C (t, M, cp), M qs (t, M, cp), and M qc (t, M, cp) through the filter component 1 16 to 
10 create the filtered signals Fj S (t, M, (p), Fj c (t, M, 9), F qs (t, M, cp), and F qc (t, M, cp), respectively. 

The analog-to-digital converter component 148 in one example converts the filtered 
signal from the filter component 116 into digital signals that represents the interference pulse 
136. The processor component 1 18 in one example samples the filtered signals F is (t, M, cp), 
Fic(t, M, cp), F qs (t, M, cp), and F qc (t, M, cp) approximately at a time t s to obtain the plurality of 
1 5 samples, as described herein. The processor component 118 employs the plurality of samples 
to calculate an in-phase term I mag and a quadrature term Q mag . For example, the processor 
component 118 calculates: 

I mag (* s . 9) = P^vy+FuK,*?) 1 = \ BT J 2(M)' cos <p\ , and 

Q mag i*. , <P) = p qc (t s ><PY+F q s(fs><P) 2 = \BT J x {m\ sin q>\ . 
20 The processor component 118 employs a modulation depth Mo such that Ji(M 0 ) is 
substantially equal to ^(Mq), for example, approximately 2.63 radians, as will be understood 
by those skilled in the art. 



Where Ji(M 0 ) is substantially equal to J2(M 0 ), I ma g(ts> <p) is proportional to |cos (p| and 
Qma g (ts> 9) is proportional to |sin <p|. The processor component 118 determines a sign of an 
in-phase term I(t s , 9) and a sign of a quadrature term Q(t s , 9). For example, 

I(ts, (p) = ± I mag (t s , 9)* and 

5 Q(ts,q>)=±Qmag(ts 5 Cp). 

The processor component 118 determines the sign of the in-phase term I(ts, 9) to be equal to 
a sign of the filtered signal F ic (ts, M, 9) or Fi s (t s , M, 9) with a largest magnitude. The 
processor component 118 determines the sign of the quadrature term Q(ts, 9) to be equal to a 
sign of the filtered signal F qc (ts, M, 9) or F qs (ts, M, 9) with a largest magnitude. As the value 

10 of to changes, the magnitude of the filtered signals Fi C (t s , M, 9), Fi s (t s , M, 9), F qc (t s , M, 9), and 
F qs (tj, M, 9) may change. The processor component 1 18 in one example determines the sign 
of the in-phase term I(ts, 9) and the quadrature term Q(t s , 9) at regular intervals to 
accommodate such a change. 

The processor component 1 18 in one example employs the in-phase term I(t;, 9) and 

15 the quadrature term Q(t s , 9) to calculate the phase angle 9. For example, the processor 
component 118 calculates the phase angle 9: 
9 = arctangent( Q / 1 ). 

The processor component 118 in one example employs the change in the phase angle 9 
between multiple instances of the interference pulses 136, 138, and 140 to determine the 
20 change in the parameters employed by the sensors 130, 132, and 134. 

Turning to FIG. 2, a plot 202 comprises a representation of an exemplary interference 
pulse 136, an impulse function h(t) for an exemplary filter component 116, and an exemplary 
filtered pulse ENV(t). The interference pulse 136 in one example is defined by: 



env(t) = exp 



108 



where T = 35 nanoseconds and To = 37 nanoseconds. The interference pulse 136 comprises a 
width of 70 nanoseconds and a rise and fall time approximately equal to one nanosecond. 
The filter component 1 1 6 in one example comprises a Gaussian low pass filter where the 
5 impulse function: 



l,* - 5643 
h(t) = exp 



tt-x ^ 
1 fc o 



where x = 10 nanoseconds. A convolution of the impulse function h(t) with the input pulse 
env(t) comprises the filtered pulse: 



ENV(t) = \h(t - v) • env{y) • dv . 



10 The filtered pulse ENV(t) comprises a time delay of 30 nanoseconds, a flat top of width 30 
nanoseconds, and a base of width 100 nanoseconds. A frequency response H(f) of the filter 
component 116 is equal to the Fourier transform of the impulse function. The Fourier 
transform is taken over a range that is sufficiently large to cover relevant areas of the impulse 
function. For example, the impulse function is approximately zero outside of a range of 0 

15 nanoseconds to 75 nanoseconds, and the frequency response is equal to: 

75 



H(f) = jh(t) • exp(-2;r li- f () dt . 



Turning to FIG. 3, a plot 302 comprises a roll-off function S(f) of the filter component 
116 of FIG. 2 where the roll-off is equal to: 



5(/) = 10.1og(i/(/)-W)). 



20 The plot 302 comprises a roll-off of over 80 dB at 100 MHz. 



Turning to FIG. 4, a plot 402 comprises the filtered signals F qs (t, cp), F qc (t, cp), F is (t, cp), 
and F ic (t, cp) with the filter component 1 16 of FIG. 2. The plot 402 comprises a phase angle cp 
equal to 0.5 radians, a time to equal to 4.2 nanoseconds, a modulation depth M equal to 2.63 
radians, and a frequency f pgc equal to 100 MHz. The filtered signals are given by a 
5 convolution of the impulse function h(t) with a product of the phase generated carrier 120, the 
sinusoidal signals I s (t), Ic(t), Q s (t), and Q c (t), and the interference pulse 136: 

F qs (t,<P) = \h(f - v) ■ S(y,<p) ■ Q,(y) ■ env(v) ■ dv , 

0 
t 

F qe {t,<P) = \Kt - v) • S(y,<p) ■ Q c iy) ■ env(v) ■ dv , 

0 
t 

F is (t,<p)= jh(t - v) ■ S(v, q>) • I s (v) • env{v) • dv , and 
o 

t 

10 F ic {t,(p)=\Kt-v)S{v,(p)-I c {v)-en V {v)dv. 

0 

The filtered signals F qs (t, q>), F qc (t, (p), F is (t, cp), and F ic (t, cp) comprise respective filtered 
pulses with a flat top of an approximate width of 30 nanoseconds and an approximate center 
at 75 nanoseconds. The processor component 118 in one example samples the filtered 
signals F qs (t, cp), F qc (t, cp), F is (t, 9), and F ic (t, q>) at the center of the pulses, for example, t s is 
15 equal to 75 nanoseconds, to promote an increase in accuracy of the samples. 

Turning to FIG. 5, a plot 502 comprises the in-phase term I(t, cp) and the quadrature 
term Q(t, cp) as functions of the phase angle cp for the filtered signals of FIG. 4. The time t s is 
equal to 75 nanoseconds, the modulation depth M is equal to 2.63 radians, and the time to is 
equal to 4.2 nanoseconds: 

20 r mag (t s9 <p) - p ic (t s ><P) 2 +Fis(t s ><P) 2 <* |cos9>| , and 

Q mag & ,9) = p qc (t 5 ,<P) 2 +F qs (t s ,(p) 2 oc |sin <p\ . 



The processor component 1 1 8 determines the sign of the in-phase term I(ts, cp) to be equal to 
a sign of the filtered signal F ic (ts, M, (p) or F is (ts ? M, (p) with a largest magnitude. For 
example, referring to FIG. 4, the processor component 118 employs the sign of the filtered 
signal F is (t s , M, (p). The processor component 1 1 8 determines the sign of the quadrature term 
5 Q(ts, cp) to be equal to a sign of the filtered signal F qc (t s , M, cp) or F qs (ts, M, q>) with a largest 
magnitude. For example, referring to FIG. 4, the processor component 118 employs the sign 
of the filtered signal F qs (t s , M, cp). 

Turning to FIG. 6, a plot 602 comprises an error A(p(t, <p) for the calculation of the 
phase angle cp of FIG. 5. The error A(p(t, q>) is equal to an output phase angle minus an input 
10 phase angle: 



A<p(t,<p) = arctan 



-<p, 



where t is equal to 75 nanoseconds. The error Aq>(t, cp) is approximately on the order of 0.1 
milliradians. 

Turning to FIG. 7, a plot 702 comprises a pulse height deviation for the filtered 
15 signals of FIG. 5. Accuracy of the calculation of the phase angle cp is dependent on a relative 
flatness of the pulses of filtered signals. A pulse height deviation AI(t, cp) is equal to the 
magnitude of the in-phase term I(t, q>) minus a magnitude of the in-phase term I(t, cp) at the 
time t s . For example, AI(t, 0.5) = 1(75, 0.5) - I(t, 0.5). A pulse height deviation AQ(t, cp) is 
equal to the magnitude of the quadrature term Q(t, q>) minus a magnitude of the quadrature 
20 term Q(t, <p) at the time t s . For example, AQ(t, 0.5) = Q(75, 0.5) - Q(t, 0.5). The pulse height 
deviations AI(t, (p) and AQ(t, cp) comprise a deviation of less than 0.1% within a flat zone of 
width 20 nanoseconds centered at t s . The pulse height deviations AI(t, cp) and AQ(t, cp) 



comprise a deviation of less than 1.0% within a flat zone of width 30 nanoseconds centered at 

The apparatus 100 in one example employs one or more computer-readable 
signal-bearing media. An exemplary computer-readable signal-bearing media for the 
5 apparatus 100 comprises the recordable data storage media 150 of the processor component 
118. For example, the computer-readable signal-bearing media for the apparatus 100 
comprises one or more of a magnetic, electrical, optical, biological, and atomic data storage 
media. In one example, the computer-readable signal-bearing media comprises a modulated 
carrier signal transmitted over a network comprising or coupled with the apparatus 100, for 
10 instance, one or more of a telephone network, a local area network ("LAN"), the internet, and 
a wireless network. 

The steps or operations described herein are just exemplary. There may be many 
variations to these steps or operations without departing from the spirit of the invention. For 
instance, the steps may be performed in a differing order, or steps may be added, deleted, or 
15 modified. 

Although exemplary implementations of the invention have been depicted and 
described in detail herein, it will be apparent to those skilled in the relevant art that various 
modifications, additions, substitutions, and the like can be made without departing from the 
spirit of the invention and these are therefore considered to be within the scope of the 
20 invention as defined in the following claims. 



